WIS/09/ 1 2- JUNE-DPPA 



Prepared for submission to JHEP 



Degenerate Rotating Black Holes, Chiral CFTs and 
Fermi Surfaces I - Analytic Results for Quasinormal 
Modes 



Micha Berkooz a Anna Frishman" Amir Zait a 

al Department of Particle Physics and Astrophysics, The Weizmann Institute of Science, Rehovot 
76100, Israel 

E-mail: micha. berkoozOweizmann. ac . il, anna. f rishman@weizmann. ac . il, 
amir . zait@weizmann .ac.il 

Abstract: In this work we discuss charged rotating black holes in AdS*, x S 5 that degen- 
erate to extremal black holes with zero entropy. These black holes have scaling properties 
between charge and angular momentum similar to those of Fermi surface operators in a 
subsector of N = 4 SYM. We add a massless uncharged scalar to the five dimensional 
supergravity theory, such that it still forms a consistent truncation of the type IIB ten 
dimensional supergravity and analyze its quasinormal modes. Separating the equation 
of motion to a radial and angular part, we proceed to solve the radial equation using a 
matched asymptotic expansion method applied to Heun's equation with two nearby singu- 
larities. We use the continued fraction method for the angular Heun equation and obtain 
numerical results for the quasinormal modes. In the case of the near-supersymmetric black 
hole we present some analytic results for the decay rates of the scalar perturbations. The 
spectrum of quasinormal modes obtained is similar to that of a chiral 1 + 1 CFT, which 
is consistent with the conjectured field-theoretic dual. In addition, some of the modes can 
be found analytically. 
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1 Introduction 

The AdS/CFT correspondence, first conjectured by Maldacena in his seminal work in 1997 
[1], provides a framework allowing us to probe the quantum nature of gravity on Anti-de 
Sitter (AdS) spacetimes via familiar quantum field theory (QFT) tools, and vice-versa. 

Black holes, being a striking feature of gravitational theories, have played a particularly 
important role in this duality and its implications. The reason is that black holes, which 
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have finite entropy, count the generic states in the theory with some fixed charges, whether 
in AdS or in flat space. Indeed, their study has shed light on the fundamental degrees of 
freedom of some gravitational background (in near horizon limits), and has the potential to 
provide more information on yet other backgrounds. If the black hole is supersymmetric, 
then under favorable circumstances one can more easily compare weak and strong coupling 
results. In cases without super symmetry, it is difficult to get an exact matching of weak 
and strong coupling computations, although in some case one can understand some general 
scaling properties. 

The most intensively studied case of the AdS/CFT correspondence is the duality be- 
tween type IIB String theory on AdS 5 x S 5 and N = 4 SU{N) Super Yang-Mills (SYM) 
gauge theory. It will also be the focus of the current paper. More precisely, we concentrate 
on rapidly rotating black holes in the U(l) 3 SUGRA consistent truncation on AdS$. Cur- 
rently, the most generic black holes known in this theory have equal charges under two of 
the f7(l)s, with all the other charges arbitrary [2] (relying on earlier work in [3—7]). We 
will be interested in the dynamical (in) stability of zero entropy degenerations of such black 
holes. Our results therefore lie at the interface of several questions: 

1. Degeneration of Black holes into singularities One of the successes of string theory is 
to resolve a large class of singularities which are ill defined in gravity. The simplest of 
these are of course the orbifold theories (for example non-SUSY singularities which 
are not accessible geometrically, cf. [8]). Still, there are more singularities which are 
not understood compared to those that are. 

The weak cosmic censorship hypothesis tells us it is difficult to create singularities dy- 
namically, in the sense that generic initial conditions will create at best a singularity 
cloaked by an horizon. One can turn this the other way around and study back- 
grounds in which there is a parametrically tiny horizon around a singularity. Such 
scenarios will hopefully capture the dynamics of the singularity while regulating at 
least some of the problems associated with studying singular backgrounds. 

The black holes we will discuss are of this nature. More precisely one can distinguish 
two cases - the first, which we will only touch upon briefly, is one where there is a 
finite entropy horizon and a singularity sits behind the horizon and approaches the 
horizon at some critical value of a parameter. 

The other, which is our main emphasis, is the case where the black hole entropy goes 
to zero. In this limit the volume factor of the horizon shrinks to zero size uniformly 
along the horizon, except at a lower dimensional locus where a singularity develops. 
Our paper deals with some general properties of quasi- normal modes in such a case. 

2. Quasinormal modes of fast rotating black holes 

Probing the linear (in) stability of a black hole background requires one to compute 
the allowed frequencies of the modes of perturbing fields. These frequencies are 
usually complex, as black holes are dissipative systems, and are thus called Quasi- 
Normal Modes (QNMs for short). Instability is observed when the imaginary part 
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causes the amplitude to grow exponentially, breaking the linear approximation. Such 
instabilities are useful when attempting to understand the phase structure of the 
black holes, and via the duality, the phase structure of the field theory at strong 
coupling. 

Computing these quasinormal modes analytically for AdS black holes is a difficult 
problem. It has been solved only for black holes which are small compared to the AdS 
radius, and are furthermore slowly rotating, cf. ref. [9]. The black holes discussed 
in this paper are fast-rotating near-extremal AdS black holes that are comparable in 
size to the AdS radius. Perturbations on such backgrounds are hard to study, even 
numerically [10] . Here we compute numerically the quasinormal modes of a massless 
scalar in the extremal limit for fast rotating black holes. For the near-SUSY black 
holes, we are able to compute analytically the lowest QNM, which is the one with 
slowest decay rate. 

The AdS/CFT correspondence relates the QNMs to the poles of the retarded Green's 
function in the dual CFT [11]. Through this relation, they can be used to study decay 
rates of excitations above a known ground state. One can therefore hope that having 
analytic expressions for QNMs will allow for precise matching to computations carried 
out on the field theory side. 

3. Fermi surfaces of fundamental fermions with the AdS/CFT correspondence 

Fermi surfaces, a robust feature of theories with fermionic degrees of freedom, are 
conjectured to have extremal black holes as their gravity duals. These constructions 
are rather constrained though. Fermi surfaces are expected to have zero entropy - 
this is true at least for the Fermi surfaces encountered in most CM systems (with the 
caveat that these are not large N systems in contrast to the setting of the AdS / CFT 
correspondence). Thus, extremal black holes, when taken to the limit in which they 
degenerate to have zero entropy, would make better candidates for Fermi surfaces. 

In this paper we propose a duality between known extremal black hole solutions in 
gauged SUGRA and Fermi surface operators living in a sector of the full J\f = 4 theory. 
Several works have suggested using Fermi surface-like operators as components of 
known 1/16 BPS charged and rotating black holes [12, 13]. It has been shown that 
while some of the properties of the Fermi surfaces are reproduced by the black holes, 
others are not. Specifically, these have N 2 worth of entropy, an indication that they 
are mixed states of many operators, and not a single Fermi surface. The zero-entropy 
degenerations that we discuss here remedy this problem. 

Using these insights, we make use of recently found black hole solutions, which we 
suggest correspond to specific operators from a restricted subset of the CFT fields. In the 
spirit of the operators discussed in ref. [13], we attempt to construct the Fermi surface 
operators while limiting ourselves to a specific sector of M = 4 SU(N) SYM. It has been 
shown (cf. ([14])) that this CFT has several sectors, which contain operators mixing only 
with operators inside the sector. The operators in these sectors, which are built from a 
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specific subset of the theory's partons, are characterized by specific relations between their 
classical scaling dimension and charges. Specifically, one of these sectors, the PSU(1, 1|2) 
sector, which we shall describe later in more detail, has interesting properties. First of all, it 
contains a small number of partons, making it easier to consider all the relevant operators. 
Second, it has an outer SU(2) symmetry rotating the two fermionic partons in the sector. 
Using this symmetry, it is possible to build operators which do not mix with any other 
operators, even within the sector. Such operators have specific scaling relations between 
the charges, which we attempt to reproduce in the gravity theory. We have identified black 
holes in the gravity theory carrying the same charges as the PSU(1, 1|2) Fermi surface 
operators. These black holes can be degenerated to have zero entropy at extremality. 
Using the conjectured duality, we can predict some of their properties from the CFT side 
by using weak-coupling computations, which we can then try and verify in strong coupling 
in gravity. 

The outline of the paper is as follows. In Section 2 we briefly summarize our results 
concerning the quasinormal modes for an uncharged massless scalar on the black hole back- 
ground. Section 3 describes in detail field theory construction of Fermi surface operators, 
the supergravity theory and the specific black hole background examined in this paper. In 
Section 4 we review the linear perturbation theory for a massless scalar, and describe such 
a scalar in the supergravity theory. In Section 5 we solve the radial equation of motion 
at the near-extremal limit using a matched asymptotic expansion method, and Section 6 
completes the analysis by solving the angular equation using an algebraic method which 
provides numerical solutions as well as some analytic results. Finally, we show some explicit 
analytical and numerical results in Section 7, and discuss them in Section 8. 

2 Summary of Results 

The black holes we are interested in are found within a consistent truncation of type IIB 
Supergravity on AdS§ x S 5 which contains the graviton, three U{1) vector fields and two 
uncharged massive scalars [2-7, 15, 16]. We examine them in the limit where the black 
holes degenerate to zero-entropy configurations, in which a ring of singularities is no longer 
cloaked by the horizon. We review the conjecture that for a specific choice of charges these 
configurations are dual to Fermi Surface operators in the PSU (1,1\2) sector of M = 4 
SU(N) SYM. This conjecture leads to predictions on the spectrum of quasinormal modes 
associated to these singularities which we attempt to verify. 

For simplicity, we shall focus here on the spectrum of scalar quasinormal modes. Us- 
ing the symmetries of the metric, the problem of finding scalar QNMs reduces to a two 
dimensional PDE eigenvalue problem. A further simplification occurs if the scalar field 
is massless, in which case the linearized equations of motion can be reduced to a pair of 
coupled ODEs, an angular and a radial one, which we shall attempt to solve. 

It should be noted at this point that in order to compare the resulting spectrum to 
the spectrum of the dual CFT, the black hole background together with the added scalar 
field should form a consistent truncation of the complete gravitational theory, i.e. type 
IIB theory. In this work we highlight the only massless scalars which can be added to the 



- 4 - 



system such that together they form a consistent truncation. Using the results in ref. [17] 
we see that coupling the black hole to a dilaton and an axion yields a consistent truncation. 
Both of these are massless scalar fields. At the linearized level, they are both minimally 
coupled to the black hole metric. 

For a generic choice of the black hole parameters, the separated equations of motion are 
both Heun equations, second order regular differential equations with four regular singular 
points. In the near-extremal limit, the distance between two of the regular singular points 
in the radial equation goes to zero. Using a matched asymptotic analysis described below, 
we are able to solve the Heun's equation with two nearby singularities. The procedure is 
then applied to the radial equation. The solution yields the separation constant A in terms 
of the complex frequency of the scalar, and we obtain 




■a rp [ 1 aX{a,m 4> ,a) 1 , 

expressed in terms of the temperature T of the black hole, which approaches in the near- 
extremal limit. Here, a is the rotation parameter around one of the independent rotation 
planes, is the azimuthal quantum number associated with the scalar field mode, and 
is the horizon's angular velocity. We have defined the parameter a = 7^^w, with 
the angular momentum and Q the U(l) charge. The AdS radius has been set to 1 
for convenience. The positive integer n is the integer number obtained from the radial 
equation. 

Next, we use the Continued Fraction Method (CFM) to solve the angular equation 
numerically for the separation constant, which together with the above expression yields 
the quasinormal modes. In addition, there are specific values of the parameters where one 
is able to extract analytical solutions for the lowest eigenvalue. The lowest eigenvalue for 
a = 2, which is close to the 1/16BPS bound, and small enough azimuthal number is 
given by 

In addition, it is possible to extract analytic results for all values of a in the large-charge 
limit, a — > 1 limit, when the black hole is rotating very rapidly. We also discuss some of the 
properties of the eigenfunctions associated with these modes, and compare them to known 
results obtained using the near-horizon techniques of Kerr/CFT. 

In order to check if these results agree with our conjecture for the dual of the black 
holes we recall that the construction of the dual Fermi surface was made in the large charge 
regime, Lea-} 1. Taking this limit in equation (2.2) we find that these quasinormal modes 
look like those of the BTZ black hole 



■( L ) =k - i4irT L (h L + n) (2.3) 
(R) = _ k _ i^T R (h R + n) . (2.4) 



It is well known that the BTZ black hole is dual to a 1 + 1 CFT and that the quasinormal 
modes above are the poles of the retarded Greens function of an operator with dimension 
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Qil, ha) in this CFT. Therefore, our result (2.2) implies that the singular degenerate black 
hole, together with the perturbation we have considered, are dual to the left sector of a 
1 + 1 CFT. Furthermore, the dimension of the operator corresponding to the perturbation 
is equal to 1, which is compatible with a free fermion bilinear. Actually, the PSU(1, 1|2) 
sector can be thought of as a 1 + 1 CFT with the coordinate (f) playing the role of x 
and thus the role of k. Thus, these findings lend support to the conjecture that the 
degeneration of a charged rotating black hole considered in this paper is dual to a Fermi 
surface construction in the PSU(1, 1|2) sector. 

3 Extremal Black Holes and Fermi Surfaces 

Several works [18-27] have suggested that extremal charged black holes are dual to Fermi 
surfaces on the CFT side. In these studies the authors compute the Green's function of 
Fermionic composite operators via a fermionic propagator in the presence of the black 
hole. They find sharp features in these functions at finite momentum and small frequency, 
which is an indication for the presence of a Fermi surface. At the same time, some of these 
features also suggest these are exotic Fermi surfaces and in particular non-Fermi liquids. 
In addition, it has been problematic to establish that these are indeed Fermi surfaces on 
the field theory side since in many cases the field theory dual is not known. For a detailed 
discussion of the subject see [28] and references therein. 

One of the subtle issues associated with this construction has to do with the entropy. 
The black holes used have finite entropy at zero temperature, whereas a conventional CM 
Fermi surface ground state would have zero entropy at zero temperature. This might be 
an effect of large-N versus finite-N, since CM systems are finite-N systems, but in any case 
it implies that the phenomenology of these black holes as Fermi surfaces would be different 
than that of known Fermi surfaces in CM systems. 

We are thus led to consider, on the gravity side, zero-entropy extremal black holes. One 
of the ways to achieve such a zero-entropy construction is to degenerate an extremal black 
hole into a zero-entropy singularity. Such an object would correspond to a system with 
zero entropy at zero temperature, at least in the large-N limit. Extremal black holes and 
their degenerations have been studied in various supergravity theories. The most studied 
solutions are in AdS 5 x S 5 and the dual M = 4 SYM [2, 3, 16, 29]. Similar constructions 
were examined in [30] for a different type of compactification, T 4 x S 1 . There, the relation 
between zero-entropy gravity configurations and Fermi surfaces in the low-energy dual CFT 
was employed in order to study the phenomenon of superradiance. 

Recently, Extremal black holes with Vanishing Horizon (EVH) have been studied in 
detail in several works [31-33] 1 , both from the gravity and field theory point of view. It has 
been shown that their near- horizon is typically a pinched AdS%, which becomes a pinched 
BTZ for near-extremal black holes. It will be interesting to examine our results in light of 
these works. 

Our goal is therefore twofold. First, we wish to identify a Fermi surface-like state in the 
CFT. Next, we would like to study its gravitational dual and verify that its extremal limit 

x We thank M.M. Sheikh Jabbari for bringing these works to our attention. 
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indeed approaches a zero-entropy configuration. We do the former in 3.1, and introduce 
our candidate black hole solution in 3.2. The rest of the paper is devoted to studying the 
quasinormal modes of a scalar in the background of the black hole and to the verification 
that they exhibit the expected behavior. 

3.1 Fermi Surface Operators and the PSU(1,1\2) Sector 

We would like to construct the Fermi surface within the context of M = 4 SU (N) SYM. 
There have been some attempts at constructing such Fermi surface operators from the 
partons of N = 4 SYM and finding black hole solutions which reproduce some of their 
properties [12, 13]. However, constructing a state using the full field content of N = 4 is 
difficult, since the theory contains several scalar fields, which usually trigger instabilities. 
Even if no instabilities exist, the scalar fields tend to mix with the fermions, potentially 
increasing the entropy. 

Fortunately, there exist several subsets, or sectors, of the theory's partons mixing 
only among themselves under the dynamics of the theory. The subsectors of N = 4 have 
been studied extensively, and have been classified in [14]. A sector which seems especially 
suited to our purposes, is the PSU(1, 1|2) sector, which contains only two fermions and 
two scalars. This sector has been analyzed in great detail in ref. [34-36]. Within the 
PSU(1, 1|2) sector, one can build Fermi surface like operators which are eigenstates of the 
dilatation operator, i.e. they do not mix with any other operators. We will consider an 
SU(4) rotated version of this sector as it appeared in ref. [14]. 

The PSU(1, 1|2) sector contains only partons with the following relations between the 
charges 

A = 2J L + Qi + Q 2 + Qs = 2J R + Q1+Q2-Q3 (3.1) 

where Ao is the classical scaling dimension, Jl and Jr are the SU(2) x SU(2) quantum 
numbers, and Qi, Q2 and Q3 are the SU (4) R charges spanning the SO(2) x SO{2) x 50(2) 
Cartan subalgebra of 50(6) = SU{4). 

The constraints (3.1) retain only four types of partons from the full set 

4 k) = vy 2i 4 fc) = vy u 

V^P^u ^=^1 (3.2) 

Here, is the covariant derivative D a a with a = a = 1. The symmetry inherited from 
the full theory is PSU(1, 1|2). In addition there is an SU(2) automorphism under which 
both $1 and (f)^ are singlets for all k, and for each k, ip^ and 'tp^ are a doublet. This 
symmetry exists only within the sector, meaning it is not inherited from the full theory. 
Following the notations of ref. [35], we use indices in order to denote the accidental SU (2) 
symmetry, defining 

v; (fc) = ^<. (3.3) 

While highly restricted, working with the full PSU(1, 1|2) sector can still be problem- 
atic. There are two scalars in the sector, which may cause instabilities and produce large 
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mixing effects. In order to restrict their mixing, we use the outer SU{2) automorphism, 
which we will refer to as pseudospin, denoted SU(2) P . Prom the charges listed above, one 
can easily see that operators which contain ip^-ip> wm generally mix with operators such 
as 0102 • However, mixing occurs only through the singlet channel of the pseudospin, 
i.e. if the fermion bilinear is taken to be in the triplet of the pseudospin, it will not mix with 
the scalar bilinear. Thus, by using fully symmetric combinations in the indices {<, >}, we 
can build Fermi surface operators out of the ^'s which do not mix with operators built 
from scalar partons. 

The simplest such construction is obtained if we build the operators from a single 
fermion In fact the collection of D^tpt, with any number of derivatives m, comprises a 
closed subsector as well, namely the fermionic SU(1, 1) sector discussed in [35] (there it is 
called the fermionic sl(2) sector). This can be shown either using the oscillator formalism 
or, in the notation we have set above, by considering the pseudospin symmetry. The single 
type of parton within this fermionic sector allows one to easily construct a Fermi surface 
operator. We review this construction explicitly below. 

This operator does not mix with any other operators in the theory, and is thus an eigen- 
state of the dilatation operator. Such an operator, containing derivatives of the fermion 
ranging from to some large K, was studied in ref. [12, 13]. Interestingly, it was found 
in ref. [13] that the dimension of such operators, both in weak and strong coupling, is the 
classical dimension with corrections of order of the inverse of the (large) charge. 

The next step is to compute the charges of these states (see eq. (3.8) for a similar com- 
putation), and look for black-hole like states with the same charges. It turns out, however, 
that they do not have finite entropy, as noted in ref. [13]. Therefore, as mentioned earlier, 
they do not make ideal candidates for studying CM Fermi surfaces. Nevertheless, both in 
weak and strong coupling the dimensions are the classical dimensions with corrections of 
order of the inverse of the (large) charge. 

One can instead use the SU{2) P to rotate the above operator to a different state in 
the same total SU{2) p multiplet. Such an operator will have both ip < and ip > and thus 
also different charges. One can then hope to find a black hole with the appropriate charges 
which has zero entropy in the zero temperature limit. We describe below the construction 
of such an operator and the corresponding black hole, which can be degenerated to zero 
entropy in the extremal limit. 

3.1.1 Explicit Construction 

The quantum numbers of and ip 1 ^ play a central role in our construction. Thus, to 
preserve the clarity of the discussion, we add a tabulation of them. 





3L 


3R 


Qi 


Q 2 


Qs 


3p 




(n + l)/2 


n/2 


1/2 


1/2 


1/2 


1/2 




n/2 


(n + 1) /2 


1/2 


1/2 


-1/2 


-1/2 



Table 1: Quantum numbers for the fermionic content of the PSU(1,1\2) subsector, j p refers to the 
angular momentum of the SU(2) p automorphism. 
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Note that each of these fields also sits in the adjoint of the gauged SU (N) so it has 
an additional adjoint index which we have so far suppressed. We will sometimes wish to 
write this index explicitly denoting the field in the following way V^q) where j = l...g and 
9 = N 2 - 1. 

We will begin by constructing the operator discussed in [13] and then present its rotated 
version, which is the one relevant to the black hole we consider below. The operator is 
given by 

2K 2K 

°i K) = n < ( i)Vi(2) • • • < (9) - n ^ a2 --- a3 €(a 1 ^> {a2) ■ ■ ■ < iag y (3-4) 

3=0 3=0 

Each term in the product is gauge invariant by itself as can be seen from the right most 
equality. To avoid cumbersome notations we will denote 

g 

Jdetm = e aia2 - a ^ (ai) * (a2) ...* K) <f = *Q) Ta (9 = dvmG). (3.5) 

a=l 

The operator described in eq. 3.4 then reads 

2K 

= Yl Jdettyi]. (3.6) 
j=o 

Proceeding, we construct the rotated Fermi surface operator with an equal number of 

i^y and ip <1 

~K-l 2K-1 

= Sym Yl Jdet i><) II Jdet 

j=0 m=K 

where Sym[ ] stands for a symmetrization of the operator with respect to the {<, >} 
indices, placing the operator in the highest SU(2) p state, with J 2 = 2KN 2 (2KN 2 + 1) and 
J z = 0. 

The operator we are considering is not a BPS operator so there is no a priori reason 
for its dimension in the strongly coupled theory to be close to the classical one. However, 
motivated by [13], we might hope that the corrections to the dimension scale as 0(-^). 
Therefore we require a large charge limit, K » 1 where we may have some control over 
the corrections. One can easily compute the charges for such an operator 

Qi = Q2 = Q = N 2 K 
Q 3 = o 

(K-\ . 2K-1 \ 

V 3=0 m=K j 



(3.7) 



K-l . 2K-1 \ 
y j=0 m=K J 

A = 2J R + Q 1 + Q 2 = N 2 (2K 2 + 2K) (3.8) 
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with Ao the classical scaling dimension. 

One could also construct similar states with higher angular momentum by leaving 
unfilled 'holes' in the operator (3.7). We are still interested in the highest outer SU{2) 
state so it should be symmetrized with respect to {<,>}. In order for the symmetrized 
state not to vanish, any ip > must appear with a different number of derivatives compared 
to any ij; < before symmetrization. For example, adding one derivative to the last operator 
in the sum, and creating the operator 



O = Sym 



A'-l 2K-2 

[ Jdet [V> <j ] ' ' Jdet [^ >m ] Jdet [^ >2K ] 

j=0 m=K 



(3.9) 



Such 'sparse' operators can only mix with other sparse states with the same number of 
partons, provided they have the same angular momentum, i.e. that the total number of 
derivatives, which are distributed among the fermionic operators, sum up to the same 
number. 

These operators satisfy the scaling relation 



a 



(Jl+Jr) 
W 



Q 

A' 2 



(3.10) 



2A- 



where for the operator (3.7) one gets a = 2 while the sparse operators have a = 2 + 

Motivated by these results, in Section 3.2 and onward, we will focus on extremal black 
holes with Q\ = Q2, Q3 = , Jl = Jr which also satisfy the scaling relations 3.8. As a 
consistency check we will verify that they indeed have E = Ao up to O (t^)- 

3.2 Description of the Black Hole Solution 

We work with a consistent truncation of type II B Supergravity on AdS§ x S 5 described in 
ref. [37]. The field content comprises of the metric, two neutral scalars and three abelian 
U(l) fields. It can also be consistently coupled to a dilaton and an axion, as we shall 
discuss later. The action of the bosonic part of the supergravity is given by 



S 



d 5 Xy 



R 



a=l 



+ E 

i=l 



Ar 2 X7 l 



+ 



/ 



fj5 J_| , \ r uvpa\ T~i jrj A k 
u 24 uv per \' 



(3.11) 



Here, l ac i s = I is the AdS radius, A 1 are the three U(l) gauge fields, and Xi are the three 
uncharged scalars, constrained by A1X2A3 = 1. These scalars may be parameterized by 



Xo 



A', 



(3.12) 



We will use black hole solutions in five dimensions that first appeared in [15], while 
following the notations of [2] which contains a generalization of these solutions. The black 
hole presented here is rotating around the two independent directions, ip and (j), with 
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angular momenta and Jj,. In terms of the field-theory angular momenta Jl and Jr 
these are 



J<p = Jl + Jr J,p = Jl — Jr- (3.13) 

This solution places no restrictions on the values of J^ and J^. The black hole is also 
charged under the three C/(l)'s. The relation between these charges and the field-theory 
charges is 

Qi = % (3.14) 

The charges J^, J^, Q\{= Q2), Q3 and E are parameterized by 5±, £3, ?n, a, 6, in the 
conventions of [2]. Our / is denoted there by 1/g. 

The near-horizon limit of these black holes at extremality was analyzed in [29] in the 
context of Kerr/CFT. The case of three equal charges was studied in greater detail [38-40]. 
Quasinormal modes when all three charges are equal were modes explicitly in [41]. 

Restricting ourselves to Q3 = and Jl — Jr = J^ = amounts to setting 6 = and 
63 = 0. The solution then reduces to the following metric 

9 x (3-15) 

Y [dt + (x 2 + 2ms*) da] \ 
(x 2 + y 2 )F 2 +V xa Xj 

where we define c\ = cosh5\,si = sinh5\ and S a = 1 — a 2 /I 2 . The functions used in the 
metric are 

X = -2m + (a 2 + x 2 ) + l~ 2 (a 2 + 2ms\ + x 2 ) (2ms 2 + x 2 ) 
Y=(a 2 - y 2 ) (1 - TV ) 



H 1 = l+ 2mS " 



(3.16) 



x 2 + y 2 



In addition, the gauge field and scalar backgrounds are 

A i = A 2 = 2ms 1 cx(dt - y 2 da) 



(x 2 + y 2 )Hi 



A z _ 2ms\y 2 dx 



(x 2 + y 2 ) ( 3 - 17 ) 



Y Y U ^/^ 

Al — A2 — ll 1 

x 3 = hI /z . 



One is often interested in the coordinate frame which is non-rotating at infinity (denoted 
by t, (ft, ip for the time and the two azimuthal coordinates, respectively), since the CFT 
quantities are defined in this frame. In the specific case we are examining, the linear 
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transformation between the different coordinate sets reduces to 

1 a , 



a = ikr T -^ (3 - 18) 

X = -ip- 

a 



Finally, using the transformation 



2 2 4 2 

x = r msi 



3 1 (3.19) 

y 



2 2 2/) 

- a cos 



in the non rotating frame, the metric explicitly takes the form of a black hole with an 
horizon of S 3 topology located at r = ro, which is the largest root of A(r) = 
The Ads/CFT correspondence suggests that the following relation holds [42]: 

Using this relation, the black hole's global charges and thermodynamic quantities are ex- 
pressed by 

vr 2ma(l + s\) 2 3 ma(l + sf) 

J (/> — o — -< * ' o 



vr m[(2(/~ 4 a 4 + S a + 1) + r 2 a 2 (S a - 2))s 2 + S a + 2] 



4G 5 £ 2 

2 2y^r 2 - fmsf (a 2 + fms 2 + r^) 2 g yAo - |msf(a 2 + \ms\ + r 2 ) 



5 = -4 — : = N' z r 6 Tt 



T 



4G5 s a s a 

r (d r A) Jrl~ fmsi 2 



47T (a 2 + r 2 , + |msi 2 ) (r 2 — |msi 



2^ 



(3.21) 



Here, ro is again the radius at which the event horizon is located. Finally, the angular 
velocity around the (ft direction at the event horizon, with respect to an asymptotically 
non-rotating AdS space, is given by 

«(i + r 2 (n, + 
n<t> = — v 9 2 1 1 ■ (3.22) 

2 1 J2. 1 2ms i 



3 



3.2.1 The extremal limit of the black hole 

In this paper we focus on the black hole very close to the extremal limit. Its temperature 
is thus very close to zero. To examine this case we need to find the location of the horizon 
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which is given by A = 0. There are two solutions for A = 0, given by X = or x 2 = 0. 
The solutions to these equations are: 



4msf 

x 2 = : r\ = — — i (3.23a) 



_1 _ a 2 _ ±£l + v / (a 2_ 1) 2 + 8m(1 + s ^ 



X = : r| = 2 !LL i i ±£ (3.23b) 



r l = - 1 -° 2 -^-^f- 1 >' 2 + 8 '"' 1 + 8 ' ) <0. (3.23c) 

Here, and for the rest of this paper, we set 1 = 1. Factors of I can be restored using 
dimensional analysis. In this notation we have: 

X = (r 2 - r 2 2 )(r 2 - r 2 ) = (x 2 + (r 2 - r 2 ))(x 2 + (r 2 - r 2 )) = (z - z 2 ){z - z x ). (3.24) 

There are two horizons r\ and r 2 and the extremal limit is obtained when r\ = r 2 in which 
case we find that the entropy is zero and the solution is not a real black hole. This is also 
the case for r\ > r 2 since x circles contract at the radius n . Only for z 2 = r 2 — r\ > we 
get a real black hole and we can use z 2 as a measure of the distance of the black hole from 
an extremal one. 

Note, that the solution has an ring of singularities at (t, x, y, (ft, ip) = (t, 0, 0, (f>, tp) which 
is hidden behind the horizon above extremality, for z 2 > 0, and becomes naked when z 2 = 0. 
Moreover, it turns out that the distance to this singularity, i.e the distance from a finite 
radial position to the horizon in the submanifold y = 0, is finite, even as the limit of z 2 — > 
is taken. Hence, there is no standard near- horizon extremal limit. 

j //y2 

It will be convenient to use a parametrization in terms of a = ^ffi^zp and z 2 . The old 
parameters in terms of a and z 2 are 

(a 2 + z 2 ) (1 + z 2 ) 2a 2 a 3 + 2az 2 

m = + — + 

I or a 

ms{ = - (3.25) 

4a + (a 2 + 1 + z 2 ) a 
zi = '—. 

a 

The BPS limit of the black hole is 

E = + .fy + Qi + Q 2 + Q 3 = J<t> + 2Qi (3.26) 

and it is satisfied for z 2 = (since the black hole must also be extremal) and a = 2. Coming 
back to the issue of superradiance, we would expect qualitatively to have superradiance 
and therefore also an instability if Vt^ > 1. Rewriting 0^ in terms of a and z 2 we get 

„ a + 2a + z 2 , 
^4, = — — . (3.27 

aa + 2 + z 2 y J 

Thus, in the extremal limit, > 1 is equivalent to a > 2 and for the BPS black hole 
Sis = 1- Therefore, for a black hole very close to being extremal, superradiant instability 
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should appear a little above a = 2. One should note that for any z 2 and a, if a = 1 then 
CiU = 1 ■ It is illuminating to rewrite the temperature and entropy of the black hole in terms 
of Z2 and a: 

s ,r2 ^^{a 2 + ^ + Z 2 ) 

t ^ 9 \ ( 3 - 28 ) 

y/z2 (4a + a + a a + 2az 2 ) 

Aair + 2-7ra (a 2 + z 2 ) 

Our main interest are rapidly rotating near extremal black holes which means we are 
considering T, a « 1 and z 2 « 1. From the expression for the entropy it is clear that the 
relative size of these quantities is quite important. In particular, in this paper we work in 
the limit ^fz~2 << S a which implies vanishing entropy. This case seems more relevant to 
a Fermi surface description. We will analyze the opposite case, for which the entropy is 
finite and large in an upcoming paper [43]. 

4 Scalar Perturbations and Quasinormal Modes 
4.1 Quasinormal Modes of Linear Perturbations 

There are several types of instabilities of black holes. The simplest one has to do with 
linearized stability about the black hole background. This will be the focus of our paper. 
The issue of thermodynamic instability and global instability to decay into other phases 
will be discussed in a separate paper [43] . 

In principle the linear analysis is straightforward, although sometimes involved - one 
writes down, and solves, the linearized equations of motion around the black hole back- 
ground for the different fields present in the (super) gravity theory. To find the energy 
spectrum of the field on this background, one uses the ansatz ty(x) = e~ tulT ^(x). In this 
convention, placing suitable boundary conditions at the horizon and at infinity implies a 
quantization of w's. The spectrum of w's and the associated modes are called QuasiNormal 
Modes (QNMs for short). For a review of QNMs in black hole physics see ref. [44]. Since 
the black hole is dissipative, most us would have an imaginary part and instability is 
associated with a positive imaginary part of u. 

Although one can write the equations and attempt to solve them numerically, this 
is usually a complicated task. It is also difficult to obtain analytical results. This is 
particularly true in AdS background where the change in asymptotic behavior introduces 
more singular points (in the technical sense of differential equations) into the equations 
of motion. Analytical results can be obtained only if the equations are separable, which 
happens only in special cases. This is also the case in which the numerical approach, 
perhaps combined with partial analytic results, is practical. 

In addition to being important quantities, characterizing different gravitational back- 
grounds, and specifically their stability, one can also hope to use QNMs in the context of 
the AdS/CFT correspondence. Roughly speaking, black holes are dual to thermal states 
in the gauge theory, and quasinormal modes would then correspond to decay times of dif- 
ferent types of excitations. However, in order to make use of the full power of the duality, 
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the perturbations studied on the gravity side must be part of a consistent reduction of 
type II B supergravity on AdS§ x S 5 . This compounds the problem since the equations of 
motion for a given set of fields in type IIB may or may not factorize, and we do not have 
the freedom of changing their couplings. 

4.2 A Consistent Truncation of Type IIB Including a Massless Scalar 

The full spectrum of the Kaluza-Klein compactification of type IIB Supergravity on S 5 
was found in refs. [45, 46]. Its truncation to the massless sector is a consistent truncation 
and results in five dimensional M = 8 gauged SO (6) Supergravity [46, 47]. The latter 
five dimensional theory consists of gravity, 15 gauge fields in the adjoint of SO(6), twelve 
2-form gauge potentials in the 6 and 6 of SO(6), 42 scalars in the 1 + 1 + 20' + 10 + 10 
representations of SO(6) (the 20' is spanned by symmetric unimodular matrices) and their 
fermionic superpartners. This theory has, in addition to the SO(6) gauge symmetry, a 
global SL(2, R) symmetry inherited from the ten type IIB dimensional Supergravity. 

A smaller consistent truncation can be constructed by restricting the above theory to 
the singlets of U(l) 3 x SL(2, R) subgroup, where U(l) 3 is the Cartan subgroup of 50(6). 
This amounts to keeping only 2 scalars out of the 20' and the three U(l) gauge fields in 
the Cartan subgroup. The resulting theory is an M = 2 supergravity with the Lagrangian 
(3.11). This theory was shown to be a consistent truncation of type IIB Supergravity in 
[37] by an explicit construction of the full nonlinear Kaluza-Klein ansatz. The black hole 
described above is a solution of this theory. 

We are interested in finding a field which can be used to perturb the black hole back- 
ground. In addition to keeping its perturbation equations simple, we require that, when 
added to the above theory, it forms a consistent truncation of type IIB Supergravity. In- 
deed, relaxing the condition for SL(2, R) invariance one gets the previous truncation with 
the addition of the axion \ and dilaton <j) (the scalars in the 1 + 1), and the following 
addition to the Lagrangian (3.11) 



As noted in ref. [17] this is also a consistent truncation of the ten dimensional theory. 
Indeed, this is the only single massless scalar enlargement which forms a consistent trun- 
cation. We cannot change the lagrangian (3.11) + (4.1), but remarkably the equations of 
motion for the axion and dilaton are separable at the linearized level for the black hole in 
question. 

4.3 Separating the Linearized Equations of Motion 

The equation of motion for a minimally coupled massless uncharged scalar in the black 
hole background is given by 




(4.1) 



-d x (xXd x <S>) + -d y (yYd y $) + 



x2 + y\2 



y 



2„,2 X 



(4.2) 
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We carry out the separation of variables using the ansatz 



$ = e -^+^<t>+™^R( x )S(y) (4.3) 

written in terms of the quantum numbers of the asymptotically non-rotating coordinates, 
r, <j), i/j. Rewriting (4.2) in terms of the coordinates appearing in the equation of motion we 
get $ = e ^x+it( am ^)+H- am ^ a2 ^)R( x )S(y). Plugging into the equation of motion 
(4.2), and separating variables, we obtain two ordinary differential equations. The radial 
equation is given by 

\d x {xXd x R(x)) a 2 m^ 2 ((x 2 + 2ms?) {-uj + am^) + am^ - a 2 uj) 2 

W) + ^ x _A ' (44) 

while the angular equation assumes the form 

\d y {yYd y S{y)) a 2 m/ (y 2 {am^ - uj) - am^ + a 2 uj) 2 

Motivated by the conjectured field-theory duals, we are interested in near extremal 
black holes which are rotating rapidly. This translates to the limits 

z 2 -»• 

a -> 1 (4.6) 

where z 2 is defined in (3.24) and a parameterizes rotation around <f>, as seen in the definition 
of the metric. The limit a — > 1 implies fast rotation, and z 2 — > implies extremality. 

In the following sections, we solve this equation in the limit y/z2 << 1 — -y, so the 
black holes being studied are arbitrarily close to extremality. In Section 5 we solve the 
radial equation analytically with the additional limit uj = m^O,^ + y^A, using matched 
asymptotic expansion. We obtain A as a function of uj and a new integer n. In Section 6 
we solve the angular equation in this limit and find the quantized values of A. Comparing 
these with the results from the radial equation yields the QNM spectrum in this limit. 



5 The Radial Equation 

We begin the analysis with the radial equation (4.4). For the solution to be normalizable, 
we impose a boundary condition at oo, R(x) — > 0. At the horizon we demand that there 
would be no outgoing waves, since we are in the classical regime. As we show in subsection 
5.1 below, the radial equation can be brought to the form of a Heun equation. For generic 
parameters it is a second order differential equation with four regular singular points. If 
one then takes the extremal limit z 2 — > 0, the equation becomes a confluent Heun equation, 
having two regular singularities and one irregular singularity. However, if at the same time 
we also take uj = m^Q.^ + y^A, the equation becomes a hyper geometric equation instead. 
In this limit, an asymptotic matching technique first developed in ref. [48], and refined in 
this work, can be employed. We describe this method in subsection 5.2, using it to solve 
the eigenvalue equation to leading order. In subsection 5.3 we apply this method to the 
radial equation, expressing the separation constant as a function of the other parameters. 
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5.1 The Radial Equation as a Heun Equation 

A Heun equation is an ordinary differential equation with four regular singular points which 
are usually taken to be located at r = 0, l,s,oo. We can permute these regular singular 
points using an SL(2, C) transformation. The Heun equation can always be brought to the 
canonical form 

r(r- s)(r-l)F"(r)+ (5.1) 
[7(r - l)(r - s) + 5r{r - s) + (a H + [3 H - 7 - 5 + l)(r - l)r]F'(r) + (a H l3 H r - q)F(r) = 0. 

The condition e = an + (3h — 7 — 5 + 1 is needed to ensure regularity of the point at oo. 

In the current form of the radial equation (4.4), it is not evident this is indeed a Heun 
equation, therefore we first perform the coordinate transformation z = x 2 . We can then 
rewrite X in terms of z and the parameters Z2, a, and a as 

X = {z- zi) (z - z 2 ) 

4a + (a 2 + 1 + z 2 ) a , , 

Z\ — (o.z) 
a 

In terms of the variable z, the radial equation indeed has 4 regular singular points, located 
at z = 0, z\, Zii oo. The outer horizon of the black hole is located at Z2- The physical region 
in our equation is thus z 2 < z < oo, which we would like to transform to the canonical 
interval [0,1]. The boundary of AdS is now located at r = 1 and the black hole horizon 
is at r = 0. To perform the transformation, we first choose the coordinate transformation 
r = fEff under which the singular points transform as: 

oo — > 1 
z 2 -> 

z 2 (5.3) 
0^ s = — <0 V J 

z\ 

Z\ — > oo 

In order to transform the equation to the desired form of (5.1) we also use the ansatz 

R(r) = (r- l)"r» (r - ^\ ^ F(r). (5.4) 

Matching to the desired form one finds that the values of the new parameters are solutions 
to quadratic equations: 

v = 0,2 

—uj (2a + a (a 2 + z 2 )) + am^ (2a + a (1 + z 2 )) 
^ = 2a^z- 2 (-z 1 + z 2 ) ( 5 - 5 ) 

am* 

Recall that the boundary conditions require the solution be normalizable, so it will be 
more convenient to choose v = 2 for the ansatz, which corresponds to R(x) ~ ^ at oo. In 
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addition, the boundary condition at the horizon requires the wave be exclusively incoming, 
i.e that the group velocity be negative. Therefore we choose 



. -U3 (2a + a (a 2 + z 2 ) ) + arris (2a + a (1 + z 2 )) 

\x — % . 

2ay/z2~(-zi + z 2 ) 

Plugging these results into the equation we find the Heun parameters: 



(5.6) 



.£2 



5 = - 1 + 2v = 3 
7=1 + 2^ 
e =1 + 2p 

an =l-i + v + p - i 
Ph =/i + v + p + i 
q =o h /3h 



Z'2 



'Zl 



Z'2 



'Zl 



"/' 



4aziy/z2~ 

a (A + uJ 2 + IQzi + 6iuy/z2 - 8z 2 ) y/~z~2~ 
2iam ( j > (3(2a + a) — iauj^/z^ + 3az2))] • 



6ia (21 + aa))uj + 6iaam^,y/—zi + a arrij ) y/z 2 + 



(5.7) 



In the extremal limit, which is of interest to us, we have z 2 — > 0, causing the points at 
z 2 and to merge. Unfortunately, they do not merge into a regular singular point, as is 
evidenced by the fact that all the Heun parameters other than 5 diverge. The equation 
then becomes a confluent Heun equation. Little is known of these equations and their 
solutions. 

If we also impose that oj = m^Q^ + y^A, the situation improves. The two points 
now merge into a regular singular point when z 2 = 0. At z 2 — > 0, we have T oc ■^fz 2 — > 0. 
The restriction to these values of oj is often used to study near-horizon quasinormal modes, 
for example in [49, 50], and for the solution of a confluent Heun equation, using matched 
asymptotic expansion, in [51-55]. 

The equation then has 3 regular singular points, and thus can always be transformed 
into the well-known Hyper geometric equation. We would like to use this fact to find the 
quasinormal modes very close to extremality, i.e. for small values of z 2 . One can easily 
check that all Heun parameters have a finite limit when z 2 — > 0. Indeed, using this scaling 
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ansatz and expanding to zeroth order in z-i one gets 



(2a + a 2 a) A 



S = 


3 


e = 


i 

L 


7 = 


1 - 






OtH = 


s( 


Ph = 







(4a + a + a 2 a) 

-io(2 + aa) 2 A + (-1 + a 2 ) J I + a 2 + ^a 2 m^ 



(2 + aa)(a + a(4 + aa)) 



(5.8) 



ia(2 + aa) 2 A + (-1 + a 2 ) yjl + a 2 + ^a 2 m, 
(2 + aa)(a + a(4 + aa)) 

-a(2 + aa)A(2i(a + a (4 + aa)) + a(2 + aa)A) + a(a + a(4 + aa))A 
9 ~ 4(a + a(4 + aa)) 2 ' 

5.2 Matched Asymptotic Expansion of the Heun Equation 

The boundary problem for the Heun equation in the limit of two nearby singularities which 
merge into a regular singularity, forming a hyper geometric equation, was studied by Lay 
and Slavyanov in ref. [48]. In this section we will carry out a general analysis of the Heun 
equation (5.1) and specialize to our case in the next section. 

We begin by following their method and deviate from it at a stage that will be made 
precise later. The essence of the method is a matched asymptotic expansion using the two 
regions r ^> \s\ and r€l where a solution in the s — > limit can be found. The small 
parameter s is used as a scale distinguishing between these two regions but also determines 
the domain of mutual validity |s| <C r ~ yjs| -C 1. 

The far region r S> \s\ 

The far region, r ^> \s\, corresponds to spatial infinity in the black hole background. In 
this region we can set s = in (5.1) to obtain the following equation 

(r - l)r 2 F"[r] + r((r - 1)(1 + a H + M + S)F'[r] + {ra H p H - q)F[r] = 0. (5.9) 

Here olh$h, etc. refer to the zeroth order contribution, in terms of s, of these parameters. 
As discussed before, the equation can be brought into a Hyper geometric equation form: 

r(l-r)G"[r] + [d - (a + b + l)r]C[r] - abG[r] = 0. (5.10) 

Using the ansatz 

F far [r]=r p G[r], (5.11) 
dividing by —r p+l and setting p to one of the values: 2 



Pi =2 ~ Ph + <5 + a/- 4 <?o + (-a/f - Ph + $) 2 

P2 = 2 1 (-aH -Ph + 8- V-^Qo + (~ q h ~ Ph + ^) 2 



(5.12) 



2 Here and everywhere else in the paper when we use a square root of a complex variable we are working 
in the following conventions for the choice of the branch cut: z — re 10 , — ty < 6 < it ; yfz = ^/re~ . Thus, 
the real part of the result is always positive. 
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will bring the equation to the desired form. For later use, notice that p\ and p 2 satisfy the 
condition 

pi + P2 = 6 - olh - Ph- (5-13) 

Choosing p = pi, which is an arbitrary choice, 3 the parameters a, b, d of the hyper geometric 
equation (5.10) read 

a =an + pi 

b=(3 H + Pi (5.14) 
d=l + a H + f3 H -5 + 2 Pl . 

There are two solutions to the Hyper geometric equation (5.10) around r = 1. Since our 
ansatz (5.4) already picks out the correct behavior near the point r = 1, we only need to 
worry about regularity of the solution and thus discard the singular solution at this point. 
We are left with the solution 

G = 2 F 1 [a,b,-d + a + b + l;l - r] (5.15) 

where 2 Fi is the Hyper geometric function. To sum up, the solution in the far region, to 
zeroth order in s, is 

F far (r)=r^ 2 F 1 [a H + p 1 ,p H + Pi,S;l-r] (5.16) 

where we have used our freedom to normalize the constant coefficient of this solution to 1 
relative to the near region solution. 

The near region r<l 

The near region, r < 1, corresponds to the near horizon domain of the black hole back- 
ground. In this region, in order to take the limit s — > correctly, we need to first make the 
transformation r = — s£ in equation (5.1). We then obtain 

- - !)(-«£ - 1)F"[£] - - - 1) - 6ti-sZ - s)+ 

(a H + (3 H - 1 -5 + - 1)(-0M + (-ohPhsZ - q)F[r] = 0. (5.17) 

Taking the limit s — > 0, the resulting zeroth order equation is 

- ^(1 + 0F"[£} - (7 + (1 + a H + p H - 6)0F'[£] - qF[Z] = 0. (5.18) 

Rewriting the equation in terms of the variable — £, we find that it takes the form of a 
Hyper geometric equation (5.10) with parameters 

a = - pi 

b = -p 2 (5.19) 
d =7. 

3 Had we chosen the opposite, p\ in the solution for the far region (5.16) would be replaced by p2- One 
can check that the resulting solution is the same solution we obtained by choosing p = pi written in a 
different way. 
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Similarly to the far region, one of the solutions of the Hypergeometric equation around 
zero should be discarded due to singularity at £ = 0. Thus the near region solution is given 
by 

FneaAO = CF (- Pl , -p 2 , 7; -£) . (5.20) 
Matching the two solutions 

The next step is to match the two solutions at the point r = ^—st with t ~ 0(1) > m 
the limit s — > 0. It is convenient to transform the Hypergeometric functions such their 
argument approaches zero in the limit s — > 0. 

Transforming the far region solution 

For the far region we will use the relation, valid when d — a — b ^ Z : 

2 F 1 (a, b, d, x) = r[a ^ 5 ~^ ]r[d] (l - rf-^F^d - a, d - b, c - a - b + 1, 1 - x)+ 

Z-lm-t MaA - d+a+b+hl - x) - (521) 

Using this relation, the far region solution becomes 
F far (t) = {V^~st) p2 r[ f 1 ~ 1 pL r[ ^i 7 2*1 {Ph + P2, an + p 2 , p 2 ~Pi + l, V^t) + 

r [a H + pi J r [p H + pij 

(v^ 1 V J^~£} T[ \ i a^i K + pi, + pi, pi - p 2 + 1, V=5t) ■ 
r [i3 H + p 2 j r [af/ + P2\ 

(5.22) 

Transforming the near region solution 

For the near region, taking r = y/—st amounts, in terms of £, to choosing £ = 77= • Thus, 
we will make use of the following formula, which is valid whenever a— b £ Z, and £ ^ (— 1, 0) 

2 F 1 ( a , b, d, -0 =r a r[^I"™ 2i?1 ( a ' a " d + lj a " b + lj t ) + 

Applying it to the solution we obtained earlier for the near region, expressed in terms of t 
and s, we get 

F near (t) =c((4A P1 r^-Hp 7 ] / _ 7 _ pi + ^ p2 _ pi + lf _ V~s 



-s/ r [7 + pi] r [-p 2 ] v i 



t Y 2 T\p 2 -p 1 ]T[ 1 ] 



—s 



Fi -7 - P2 + 1, -P2, Pi - P2 + 1, -■ V- ■ ( 5 - 24 ) 



r [7 + P2] r [-pi] v i 

Now we have everything set to expand both solutions in powers of s and compare the 
matching powers of t. We will use the fact that, unless — d € IN, 2 F\ [a, b,d,s] = 1 + O(s) 
for s — > to write the two equations resulting from the matching between (5.24) and (5.22). 
Equating the coefficients of the powers t pi in both equations gives 

r= (srr[p 2 -p 1 ]T[5]T[ 7 + p 1 ]T[-p 2 ] 

T[ Pl -p 2 ]T[~(\T[5-a H - Pl \T[5-p H - Pl Y l " > 
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Using this result and comparing the powers t P2 we obtain (here we deviate from the results 
in [48]) 

(r [pi - p 2 \) 2 r [7 + p 2 ] r [-pi] r [fa + p 2 ] r [ag + p 2 ] = , r—,2 P1 -2 P2 

V[a H + p x ]V[p H + p x ]{T\p 2 -p l ]yT[ 1 + p 1 ]T[-p 2 ] [V > • [ - ' 

Assuming Re(p\ — pi) > 0, the RHS of (5.26) goes to zero as s — > 0, so the denominator of 
the LHS must diverge in when taking the limit. This means one of the arguments of the 
gamma functions in the denominator should be equal to a negative integer which gives a 
quantization condition. 

Up until now the discussion was completely general. In the next section we will focus 
on the radial equation, expressing the quantization condition explicitly. 

5.3 Applying the Asymptotic Matching to the Radial Equation 

In order to extract the correct quantization condition for A from equation (5.26) we must 
check several cases, some of which do not lead to a valid solution. As we have shown, the 
matching procedure implies that one of the following should be integer: 



<*H + Pi = I 3 + X 1 1 " - + fl(4 + aa) + (2 - aa){a + fl(4 + aa)) I 




aX (-1 + a 2 ) ^/1 + a^ + f a^m 

I3h + Pi = ^ I 3 +l/ 1 - a + a(4 + aa) (2 + aa)(a + a(4 + aa)) 1 (5 ' 27) 



aX (-1 + a 2 ) ^l + a2 + ^a 2 m/ 



ia(2 + aa)A I aX \ ,„ 

7 + pi = - 1 y —- J — + Jl 5.29 

a + a (4 + aa) V a + a (4 + aa) J v ; 



ia(2 + aa)A / aA \ 

a + a(4 + aa) V a + o(4 + aa) / v ; 



aA 

p 2 -pl = -<l/l ; 7— r- 5.31 

a + a(4 + aa) 

The expressions possibly consistent with this condition are 4 (5.28), (5.29) and (5.31). In 
Appendix A we show the expression (5.31) doesn't lead to a valid quantization condition. 

Quantization condition for small 

We are interested in rapidly rotating black holes, i.e 1 — a 2 << 1. Thus for finite fixed 
the third term in (5.28), which is the only negative term with a negative real part in the 
expression, is very small compared to the other terms making (5.28) positive. Therefore, 
the only possibility for small is 7 + p\ = — n which gives the following quantization 
condition for A: 

-4n(l + n)(a + a(4 + aa)) + 2m(l + 2n)(2 + aa)A a 2 (2 + aa) 2 A 2 

A — 1 -. t- rr, (5.32) 

a a(a + a(4 + aa)) 



4 Since a < 1, the expression in (5.27) is manifestly positive, also it is easy to see that 7 + pi = —p2, so 
both expression (5.29) and expression (5.30) will provide the same quantization condition. 
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or, more compactly 



A= ° + 4 ° + a2 ° W °t , -l-,(l + 2„) . (5.33) 
2a + a z a \ V a + 4a + a z a / 

Taking into account Re {p\ — p^) > 0, which is necessary for the validity of the calculation, 
we get the additional restriction 

1 . „ 1 „ / ia(2 + aa)A \ 

Re(-n - 7 - - - ^ + 5)) = - Re -1 - 2n+ — K —— — J — > 0. (5.34) 



2 V 2 V a + a(4 + aa 

This in turn implies 

, 4a + a + a 2 a , . , „ . „ . 

ImA^< - 5 2n + l = -2irT 2n + 1 5.35 

2a + a^a 

where T is the temperature of the black hole. In particular, this means Im A < 0, so the 
chosen quantization condition gives only stable solutions. 

Thus we obtain that A is always purely imaginary which implies 

Reu = m^ft^ (5.36) 

for any parameters and small enough m^. Finally, using the restriction (5.35), we can 
eliminate the solution with a plus sign in front of the square root in 5.33. 

Quantization condition for large 

Turning to > taking an + pi = —n we obtain an expression for A which is 

independent of A. This means that at leading order the angular equation won't include 
A, and one needs to compute the next-to-leading order in yjz^ to extract a quantization 
condition for it. 

As in the previous case we apply the condition Re {p\ — P2) > to get 



.„_ l ( - aH - fa+ ,) * n _^^ + V 1 Y + ;^ >0(5.37) 
2 2 /i + a 2_|_4a 4 + 2aa 



and the restriction on m^, coming from the demand that the solution be compatible with 
n > 0, is 

3(2 + qq) yJa + a{A + aa) 
m <t> > (l- a 2) a 3/2 ■ (5 - 38j 

This is the quantitative interpretation of the requirement for to be large in order to 
ensure the validity of the quantization condition. D is a positive quantity, and so no 
instability is caused by the quasinormal modes obtained in this limit. 

To sum up, there are two valid quantization conditions one of which always gives stable 
solutions and is valid for any m^. The other condition is valid only for large enough m^s 
and one must employ the next to leading order to find the QNMs. Fortunately, in the 
limit of interest, i.e. of rapid rotation a — > 1, equation (5.38) implies that these modes are 
pushed to infinite for finite small enough Z2, and they disappear from the spectrum. 
We will therefore ignore them in this paper, which analyzes low energy excitations of the 
black hole, and we will revisit them in a later paper [43], which analyzes high frequency 
modes in a slightly different limit using numerical tools . 
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6 The Angular Equation 



The angular equation is given by 



la, te Y3 s Sfa» - ^ («-*-")-«-* + -A-) + + A \ m _ . (tl) 

As we discuss below, it can also be brought to the canonical form of a Heun equation, 
shown in (5.1). In subsection 6.1 we show the explicit transformation. However, there 
does not appear to be a limit, consistent with the Z2 — > limit we took in the radial 
equation, in which it simplifies and allows for a general analytic solution. We therefore use 
the Continued Fraction Method (CFM), described in subsection 6.2, to solve the eigenvalue 
problem and find the quantization condition. In subsection 6.3 we apply the method to 
the angular equation, and obtain numerical solutions for eigenvalues, as well a an analytic 
result in the case a = 2. 

6.1 The Angular Equation as a Heun Equation 

In terms of the variable z = y 2 the equation is also a linear ODE with four regular singular 
points: z = 0, a 2 , 1, oo. The physical region of interest is < z < a 2 , which we would like 
to transform it to the interval [0,1]. Thus, we use the transformation r = under which 
the singular points transform as: 

0^0 



a 2 ^l 



72 > L ( 6 - 2 ) 



a 



(6.3) 



Next, we plug in the ansatz 



1 



x 



S(r) = r K (r - If [r - H{r) (6.4) 



a 2 

and using any of the following choices 

» = ±f ^ = ±^ X = ±| (6.5) 

the equation takes the form (5.1). In order to have regularity in the interval y £ [0, a] 
the preferred choice in our case is r\ = -^p- and k = i^jp. The parameters of the Heun 
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equation, as denned in equation (5.1) 



1 

s = — 

or 

1 n n 
a w =k + jj + x = - (|m^| + \m^\ - u) 

0H =2 + k + i] + x = 2 (4 + |»7ty| + |m^| - w) 
J =1 + 2r? = 1 + |m^| 
7 =1 + 2k = 1 + | m,/, | 
e =l + 2x = l-a; 
_-2a 2 w - 2 ((-1 + a 2 w) |m</,| - \m^\ (1 + |m^|)) 



(6.6) 



4a 2 

^2 (2 (m^l + mj) + ^2 ( A + (~ aw + m <^) 2 + "4 

Both choices of % should produce the same results. It is more convenient however to choose 
X = — % in order to maintain the explicit dependence in the parameters on uj — for 
771^ > 0. 5 

6.2 The Continued Fraction Method for the Heun Equation 

There are two local solutions to Heun's equation (5.1) in the neighborhood of the singular 
point r = and they take the following form 6 

oo 

(6-7) 



k=0 

The relevant boundary conditions for our problem is that the function approaches a con- 
stant at r = and r = 1, as we have stripped off the angular dependence in the ansatz 
(6.4). Recall that s > 1 in our case. Plugging the series solution around r = 0, denoted 
previously by H®(r), back into the Heun equation results in a 3-term recurrence relation: 

/fcO fc +i + guak + h k a,k-i = 0, (6.8) 

with the recursion coefficients given by 

fk =a(fc + l)(fc + 7 ) 

g k = - k[s (k + 7 + 6 - 1) + k + 7 + e - 1] - q (6.9) 
h k ={k + a H - \){k + p H ~ I)- 



5 The alternative choice can be reproduced by reinstating I, as appropriate by dimensional analysis, and 
making the transformation to I — ¥ — I in these equations. 

6 The same is true for the point r = 1 with the replacement of 7 by S and of r by 1 — r. 
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The initial conditions for this series are 

o_i = 

(6.10) 

a = 1. 

Such three term recurrence relations, in the context of solutions to differential equations, 
were studied thoroughly by Leaver [56, 57] and we will follow his methods to obtain a 
numerical solution. 

Within this method, regularity at r = 1 restricts the coefficients in (6.9) which in turn 
determines the quantization of A via the dependence of q on A. 

Regularity at the point r = 1 

Indeed, for the series solution to be a genuine solution it must also satisfy the boundary 

oo 

conditions, and in particular be regular at r = 1. This implies the convergence of ^ a^. 

k=0 

To check whether the series converges we examine the asymptotic behavior of \r^\ ~ 
Denoting 

t = lim rfc = lim rk+i 
9 k 

a = lim — 

k— >oo Jfc 

b = lim ^ (6.11) 

k— >oo Jk 

we rewrite equation (6.8) in the form: 

t 2 + at + b = 0. (6.12) 

The two solutions to this equation give us two types of asymptotic behaviors for r&. 

Actually, any three term recurrence relation has two independent solution sequences 
(without stating the initial conditions) Af. , B/. . These frequently have the property that 
lim 4 41 = in which case is called the minimal solution. If \ti\ > fa] then the solution 

k— >oo fc 

with the asymptotic behavior t2 is the minimal one. There is a unique choice of initial 
conditions corresponding to the minimal solution while all other choices inevitably produce 
the dominant one. As we explain below, the correct solution to the problem is the minimal 
one. Therefore, since the initial conditions are set and are independent of any parameter, 
as seen in (6.10), we will need to impose a restriction on the recurrence coefficients (6.9) 
such that the resulting series solution is minimal given these initial conditions. This will 
provide the required quantization condition for A. 

Applying the analysis above to the Heun equation we have 

h =i 

1 (6-13) 

t 2 =-■ 

s 

Recalling our convention, s > 1 we see that the minimal solution will always give a con- 
vergent series at r = 1. 
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However, the dominant solution might also give a regular solution at r = 1. To check 
the convergence of the dominant solution we retain the ^ terms in a and b to get 

-2 + 5 

h = 1 + — (6.14) 

n 

Thus, depending on the value of 5, the dominant solution might give a convergent series at 
r = 1. However, regularity alone is not always enough and, as described in [58], using the 
full boundary condition at r = 1 it can be proven that the minimal solution is the correct 
solution of the problem for any 5. 

Taking into account the boundary condition at r = 1 

For completeness we repeat the argument of [58] here. Returning to the boundary condition 
at r = 1 we recall we have assumed that the corresponding power law behavior was stripped 
off, so for r = 1 the function should simply go to one. To find the restrictions stemming 
from this demand we look at the two local solutions around the point r = 1, they have 
the properties H\(r) ~ 1 + • • • and H^ir) ~ (1 — r) 1_<5 + • • •. Any differential equation 
of second order has only two linearly independent solutions. In the interval [1, s] the local 
solution around zero, which we considered thus far, should be a linear combination of the 
two local solutions around r = 1. Moreover, to satisfy the boundary condition it must be 
exactly equal to H\(r), which implies it should be absolutely convergent up to r = s and 
not only up to r = 1 (and also that it is equal to what is called in the literature a Heun 
function). Examining equations (6.13) it is clear this requirement is satisfied only for the 
minimal solution. Thus, we have shown that for a perturbation of a black hole described 
by the Heun equation, in the context of the CFM method, the correct solution in the whole 
range of the parameters is the minimal one. 

Obtaining the quantization condition 

Returning to the issue of a quantization condition, we quote a known result for three term 
recurrence relations. The following equation holds for the minimal solution of the recursion 
relation (6.8), and in particular the continued fraction shown below is convergent [59] 

Qfc+i _ ~ hk+i h k+ i fk+ih k+2 /fc+2^fc+3 (6 15) 

a k fk+i ■ h k+2 9k+i~ 9k+2~ 9k+3~ 

fk+2 • hk+3 

9k+2 

9k+3 - ■■■ 

Using this equation for k = together with the equation 2i = — 22. (arising from the initial 
conditions (6.10) and the recurrence relation) we get the condition: 

n /o • hi fl-h 2 f2'h 3 /aia\ 

= 5o ••• (6.16) 

9i- 92 - 93- 



7 In the example of the angular equation we will treat below, regularity is enough to uniquely choose the 
minimal solution. However we would like to show that the minimal solution is the unique solution also in 
the more general setting. 
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Truncating the continued fraction at some k we thus get a polynomial equation for A which 
can be solved numerically. Since the continued fraction is convergent, some of the As found 
by solving this equation correspond to true eigenvalues, with the smallest ones found first. 
Increasing the cutoff causes more QNMs to be found, and also results in better accuracy 
for the ones already found. 

6.3 Applying the Continued Fraction Method to the Angular Equation 

Before using the CFM to solve the angular equation we first need to take the same limit we 
took when dealing with the radial equation. In particular, we replace uj = m^Q^ + y/z^A 
and m^h = into the relations (6.6) and expand to zeroth order in zi- The parameter A 
enters these relations only through A in this order. Thus, one can proceed to solve (6.16) 
numerically for A, getting quantized values for it. Then, rewriting (5.33) to express A 
in terms of A, one can use each of the values found for A to extract another quantized 
spectrum for A. 

When carrying out this procedure, there is one special case. Setting either a = 2 or 
a = 1, after taking the limit discussed above, implies = 1 + O(^)- This in turn means 
that co = rrifj) to zeroth order in Z2 and one obtains in (6.6) 

cm = 

(6.17) 

Ph = 2. 

Note that the above relations hold also if one sets = or takes a — >■ 1. In the recurrence 
relation (6.8), hf. now takes the form 

h k = (-l + k)(l + k). (6.18) 

In particular, we have h± = . 

Plugging hi = into the recurrence relation (6.8) we find that the equation (6.15) is 
no longer valid for k = 0. Instead one obtains an initial condition for k = 1, which is 

^ = ^ (6.19) 

This condition can be used to derive a quantization condition in the same fashion it was 
derived for k = 0. The new quantization condition is 

= gi -fLJ».f*J*... (6.20) 

92- 93~ 

However, a much more interesting consequence of the relation hi = is that, if we 
assume 

90 = (6.21) 

and use the recurrence relation (6.8) at k = with the initial condition a_i = 0, we get 
that a\ = 0. Then, for k = 2, the recurrence relation implies that a,2 = which means that 
all coefficients of the series are zero except for ao which is equal to 1. Thus, the assumption 
(6.21) leads to the solution 

H(r) = 1 (6.22) 
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which meets the requirements of the boundary conditions and thus is a legitimate and 
exact solution of the problem. Moreover, (6.21) supplies the corresponding equation for 
the eigenvalue A, which can be solved analytically. 

A simpler way to see H(r) = 1 solves the equation is to examine (5.1) under the 
condition an = 0. It is evident that a constant solves this equation if both q = and the 
boundary conditions are H(0) = H(l) = const. Indeed, the condition q = is equivalent 
to (6.21) and in our case the boundary conditions are H(0) = H(l) = 1. As mentioned 
above, this also occurs when a — > 1 or = 0. 

Indeed, we obtain the value of the separation constant 

A = -J7ty(l - o)[2(l + a) + j7ty(l - a)]. (6.23) 
Plugging A into equation (5.33) 8 we find 

A = ■ : - 2 ( 1 + a )( 1 + ") + (- 1 + °W (6 24) 

a 

7 Results 

We have analyzed the equation of motion of a massless uncharged scalar in the background 
of a charged rotating black hole. The equation of motion for such a scalar separates into an 
angular and a radial part, and we have denoted the separation constant by A. The radial 
equation can be solved in the extremal limit, Z2 —> 0, along with u = m^l^ + ^/z^A and 
= 0. As noted in [40] for similar black holes, in the above limit the angular equation 
does not depend on A. The eigenvalue problem for A is thus decoupled, and one can solve 
it independently, subsequently using the result to solve the eigenvalue problem for A in 
the radial equation. 

The Radial Equation 

Solving this equation using the matched asymptotic expansion yields the following relation 
between A, and hence u, and A 



11 a\(a, a, 7tu) 1 \ .mm/ 

u = m^&fy - iAttT -4/1 2 + - + n\ = m^l^ - iAitT (D (a, a, m^j + n) 

\ — y ca ~\ qcl ~~\~ a cy. — ■ j 

(7.1) 

expressed in terms of the temperature T of the black hole, which is used here as an inde- 
pendent variable instead of zi- In this limit T = \J~z%(a + l)/(27ro). The deviation of u 
from m^Qfi must be purely imaginary in order to ensure validity of the solution. Apart 
from its dependence on m<k, this is the same spectrum as for an operator of dimension D 
in a finite temperature chiral (1 + 1) CFT. 



3 As we discuss in the paragraph below (5.33), the correct choice of sign in (5.33) is a minus. 
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The Angular Equation 

The angular equation can be used to solve for A in terms of and a. If one assumes 
either a = 2, a — > 1 or = it is possible to obtain some eigenvalues and eigenf unc- 
tions analytically, while for a generic choice of these parameters we must calculate A only 
numerically. 

Numerical Results for the Quasinormal Modes 

We now list some numerical results obtained using the continued fraction method. The 
computations have been performed using Mathematica, which allowed for high precision 
numerics. The parameters are chosen to be a = 2, = 5 and a = 0.93. The choice of 
a allows us to compare the lowest mode with analytical results. The results displayed are 
for D, with uj = — 4miT(D + n), and they are valid for T « 1 — a 2 . 



A 


-1.47 


-11.17 


-24.86 


-42.32 


-63.35 


D 


1.09 


1.50 


1.89 


2.26 


2.62 



Table 2: Computed values of the separation constant A and the offset D, computed using the CFM 
for a = 2, TO0 = 5 and a = 0.93. 

In Figure 1 we plot the first quasinormal mode for n = and a = 0.93 for several as 
and 77US. 
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1.15 



1.10 



1 .05 

















m<s=9 



1.4 



1.6 



2.0 



2.2 



2.4 



2.6 



Figure 1: The parameter D for the first quasinormal mode with n = and a = 0.93. The grey 
line marks a = 2 for which we obtained analytical results as well. 



As was previously mentioned, the singularity remains at a finite distance even when 
the extremal limit is taken, when one approaches the horizon radially along y = 0. One 
may expect that in such a case there would be no decoupled near- horizon limit, which 
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occurs when the horizon is pushed to infinite proper distance. Performing the analysis we 
find that as z<i tends to zero the number of modes grows, and correspondingly they become 
denser. This is easy to see because the spectrum of A, and therefore A is independent of z% 
in this limit, and as Z2 — > it condenses around m^Q.^. This is an indication that modes 
on this background have a decoupling limit, similar to that of an AdS throat. 

Analytic Results for a = 2 

• For the case of a = 2 we obtained an analytic expression for one of the eigenvalues 
A, given by 

A = -77^(1 - a)[2(l + a) + m^(l - a)]. (7.2) 
The analytic expression for the eigenfrequencies in this family of modes is given by 



uj =m, 



+ V*2 (-(2 + 2n) U + + ™<t> U- ^Y) 



• Using the numerical procedure, we checked several sets of parameters in which the 
analytic solution exists. Sorted by their absolute value, the mode obtained analyt- 
ically is the lowest of the bunch. Thus, in terms of uj, we get a family of modes 
parameterized by n, with n = corresponding to the mode with slowest decay. 

• The analytic expression for the angular part of the eigenfunction of the lowest stable 
mode is given by 

m d> m 4> mj, -mj, 

S(y) = (y 2 -a 2 ) 2 (y 2 - l) 2 = {a 2 sin 2 6) — {I - a 2 cos 2 6)— . (7.4) 
where the change of variables from y to 9 is y 2 = a 2 cos® . 
Results in the Large Charge Limit - a — >■ 1 

As previously mentioned, one can also obtain some analytical results upon taking the 
additional limit a — > 1 or setting = 0. It is easily seen that in this case, A = is a valid 
eigenvalue with S(y) = const, and consequently, for this mode D = 1. 

In fact, when a — > 1, the angular equation does not depend on or a. The de- 
pendence on a enters only through $7^, which approaches 1 in this limit, and then 
enters the equation solely through the combination (1 — a)m^,. Thus, we obtain a universal 
spectrum, to leading order in 1 — a, given by 

cj = nifj, — i4irT (D(a) + n) , (7.5) 

with the dependence of D on a obtained using (7.1) 



„ , ^ 1 / a\ u 1 

D(a) = -\1 ^t + -. 7.6 

v ' 2 V 2a + 4 2 y 1 
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Here, A n is the universal angular spectrum in the a — > 1 limit. We have computed its first 
lowest-lying values, which are given in the following table. 








-10.12 


-25.44 


-44.95 


-68.17 



Table 3: Numerical values for the lowest modes in the universal spectrum X u . 

Some General Features of the Radial Part of the Wave Functions 

An expression for the radial part of the wave function, as a function of the separation 
constant, was obtained in Section 5.2 for two different regions of the r segment and a 
matching between them was performed. The exact expressions for the function in the 
two different regions when a = 2, and an analytic expression is available, are written in 
Appendix B , here we only describe the most prominent features of these solutions. 

First, we consider the region stretching from infinity to a distance of — s = — — from 

2 

the horizon. In terms of r = ■> , f, v > , this is the region J — s < r < 1. We find that in 
this region the function is a polynomial in r and is independent of n. Thus, far from the 
horizon the wave function &(r,d) = e -^~ m ^) R(r)S(6) depends on n only through its 
time dependence. 

It is important to mention this is true not only for the lowest eigenvalue for A, which 
we found analytically, but also for the rest of the modes. This is a consequence of the 
fact that A is independent of n, thus if one changes A it cannot affect the n dependence 
of the eigenfunction. Furthermore, this statement can be verified also for a ^ 2 since 
the n dependence can be determined from the expression for the radial function without 
specifying A. Indeed, we find it is independent of n. The precise expression in terms of A is 
given in Appendix B. Thus, the eigenfunctions away from the horizon are all independent 
of n, and couple to all the near-horizon modes. 

In the near horizon region, < r < y 7 — s, the most prominent feature of the solution is 

£) n 

that when one approaches the horizon, r = 0, it takes the form R(r) ~ r _ n — \-0( — ) where 
D is defined in 7.1. One implication of this result is that the solution diverges for n > as 
one gets closer to the horizon. This is a general statement for decaying quasinormal modes 
in the background of a black hole. One way to see this is to examine the behavior near the 
horizon. For a scalar field one usually gets 

r ~ (7.7) 

where A and B are independent of uj and A > 0. The relation is written using coordinates 
in which the horizon is at r = 0. The boundary conditions at the horizon require a purely 
ingoing wave which leaves us with the minus sign in (7.7) (cf. [9]) . Thus for a decaying 
mode, i.e Imw < 0, since uj is the only complex quantity in the equation, r appears with a 
negative power in $ near the horizon. Indeed this phenomena appears in other works on 
quasinormal modes of black holes, such as [56] and [9]. 



(J> ^ r ±A{iw+B) e -iwT 
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8 Summary and Discussion 



This work analyzes the equation of motion for a massless minimally coupled scalar in 
the background of a charged rapidly rotating black hole in the limit of T — > 0, which 
corresponds to an extremal zero-entropy degeneration of a black hole. 

We began by examining Fermi surface operators in the PSU(1, 1|2) sector of N = 4 
SU(N) SYM theory. These operators have been constructed explicitly out of the fermionic 
operators of this sector and, due to its extended symmetry, are eigenstates of the dilatation 
operator. They are also known to have, in the large charge limit, small corrections to the 
conformal dimension. One may hope that this would allow matching them to gravity 
configurations of similar properties. 

A conjecture has then been made regarding the gravity solutions dual to these black 
holes. By studying black holes which are solutions of a consistent truncation of type IIB 
SUGRA on S" 5 , the dual configurations can be identified using scaling laws between the 
angular momentum and charge of the Fermi surface operators. In order to verify the 
correspondence, we compute the quasinormal modes of an massless scalar field. For the 
result to have a meaning in the dual CFT, the inclusion of the scalar must also reside in 
a consistent truncation of the full theory. The addition dilaton (or axion) forms such a 
consistent truncation. 

The limit T — > corresponds to a singular background, as the ring of singularities, 
located along one of the independent rotation planes, reaches the horizon. The area of the 
horizon also shrinks to zero in that limit, forming a zero entropy degeneration. Nonetheless, 
it is possible to solve the equations of motion semi-analytically for all quasinormal modes 
and for all black holes in the family we identified. Even more surprising is the fact that 
one can find a full analytic expression for the quasinormal modes of the lowest lying stable 
mode of black holes with a = 2 and T — > 0, with a = j^j^jji- Such black holes are 
arbitrarily close to the supersymmetric black hole located at T = 0. 

The Emergence of a 1 + 1 CFT and its Relation to the Conjectured Fermi 



• It is interesting to compare the leading order of the large charge limit of the quasi- 
normal modes we obtained to the quasinormal modes of a BTZ black hole. In the 
large charge limit, corresponding to a — > 1 on the gravity side, the dependence of A 
on rricf, and a vanishes to first order in 1 — a. Then, the spectrum is 



For comparison, the expressions obtained in [11] for the quasinormal modes of a BTZ 
black hole are 



Surface 




(8.1) 
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Indeed, the quasinormal modes we obtained look exactly like those of the left sector of 
the BTZ black hole, when replacing the fourier variable corresponding to momentum, 
k, with that corresponding to the angular momentum in our case, m^. 

• Recall that, as shown in ref. [11], the quasinormal modes of perturbations around a 
black hole background correspond to the poles of the retarded Green's function on 
the CFT side. In the case of the BTZ black hole the dual CFT is a 1 + 1 field theory. 
Thus, the structure of the quasinormal modes we found implies that our extremal 
black hole, in the large charge limit, matches the left sector of a 1 + 1 CFT. 

• As discussed in the previous section, in this limit, as well as when = 0, there is 
always a zero eigenvalue for A, and the spectrum corresponding to this eigenvalue is 

u = - i4vrT(l + n) . (8.4) 

Thus, the lowest excitation is an operator with dimension hi = 1 in a 1 + 1 CFT. 

Given this last observation, we would like to examine the compatibility of the above 
results with our conjectured dual Fermi surface operators. We have suggested that this 
family of black holes corresponds to the PSU(1, 1|2) sector which is limited, in terms of 
derivatives, to D^. This means that fields inside the sector are constrained to rotating 
with a specific chirality in the eft direction and effectively live in a chiral sector of a 1 + 1 
CFT, with the x direction in the standard CFT replaced by <j>. Note that the perturbation 
we considered, the massless uncharged scalar <I>, is inside the sector, as it has = 0, all 
charges equal to zero, and oj = (see relation 3.1 and 3.13). Thus, it is quite natural 
that we find poles corresponding to an operator living in the left sector of a chiral 1 + 1 
CFT, if one replaces the symmetry for translations along x by symmetry for translations 
along 4>. 

Furthermore, taking the large charge limit, we obtain the spectrum of a free fermion 
bilinear, D{2) = 1, in such a theory. Recall that $ in AdS 5 x S 5 /Af = 4 SYM is dual to 
the lagrangian of the CFT. When restricted to the PSU(1, 1|2) sector, only the fermion 
bilinears remain. Hence, it is natural for the dilaton to excite a fermion bilinear in this 
1 + 1 chiral field theory. The decay rate of these excitations is very small since it scales 
with the temperature T. Therefore these excitations are quite long lived. 

In the context of the Fermi surface operators, we are interested in values of a of the 
form a = 2 + o(g) (see the discussion in the end of section 3.1.1). Curiously, when a = 2, 
we find analytically the lowest mode for every value of a and m^. 

Our results lend support to the conjecture that this degeneration of a charged rotating 
black hole is described by a chiral 1 + 1 CFT, which is a Fermi surface of the microscopic 
fermion. To verify this, the next step would be to compute the poles in fermionic Green's 
function, as well as those of additional, charged, scalar fields. Not only should the results 
agree with those of a chiral 1 + 1 CFT, but the coupling of these bulk fields to the black 
hole are unambiguously prescribed by our model in Section 3. We hope to return to this 
in future work. 
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At first glance, our results are reminiscent of the results obtained using Kerr/CFT 
[60] and its extensions to AdS, higher dimensions and gauged supergravities, e.g. [61, 62]. 
For a review, see [49]. In fact, Kerr/CFT has been applied to a broad family of extremal 
backgrounds in gauged supergravities [29]. Our black hole is obtained as a limit of these 
backgrounds, but is not covered by these works as the extremal near horizon transformation 
used there becomes singular in our case. More specifically, our black hole has zero entropy 
at extremality, and thus it is unclear how to apply the Kerr/CFT duality in this case. As 
mentioned, the more recent work on EVH may be more appropriate [32, 33]. 

It would be interesting to examine in more detail the near horizon structure of this black 
hole when it becomes singular, which we have not discussed in this work. Understanding 
the singularity structure at y = and zero entropy degeneration at extremality can help 
shed light on the dual low energy CFT. 

Finally, The fact that we have obtained solutions in the form of Hypergeometric func- 
tions suggests the presence of some SL(2, M) symmetry in the near extremal limit, at least 
for the massless scalar in this background. It may be useful to understand how the SL(2, M) 
symmetry of the obtained solutions manifests itself in this scenario, compared to previous 
results, such as [40, 63, 64], as well as its possible connection to the AdS^ factor present in 
the near horizon limit of the black hole away from the singularity. 
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A Elimination of additional possible quantization conditions for A 



First, examining the case p\ = p2 we discover that the equality in equation (5.26) is satisfied 
without further restrictions. This implies the relation p\ = p2 is exact, without any non 
analytic corrections. However, if we wish to consider this case we can't use the results 
achieved thus far since the linear transformations we used in equations (5.21) and (5.23) 
are no longer correct and there are other relations specific for this case. Repeating the 
process with the changed relations amounts to an equation including Zn(|s|) on the one 
side and a sum of Digamma functions on the other. Since Digamma functions can only 
have simple poles, this equation can never be satisfied and therefore we must discard this 
solution. 

Thus, the remaining case is Re(/9i — p<i) > 0, for which we still haven't exhausted all 
the possibilities for equality in equations (5.27)-(5.31). Specifically, we need to check what 
happens for p2 — p\ = — n + O(s). This corresponds to the special case we mentioned 
above, where —d E IN for the Hyper geometric functions both in the far and near regions 
and so should be treated with more care. In particular, before taking the limit s — > in 
equations (5.24), (5.22), we should consider that, provided either a or b are not equal to 
some negative integer —m, satisfying m < n, the following holds 

^) = (aWi(iW V , [fl + k + 1 , b + k + 1]k + 2]z] . (A.!) 

We can now make use of p2 — pi + 1 = — k, i.e n = k + 1, along with the formula T[— n + 1] = 
— nT[— n] and divide equations (5.24) and (5.22) by T[p2 — pi]- This allows us to take the 
limit s — > 0. For the zeroth order of equation (5.24) we get 

F{ t) = ^ [Pl PA p7] i (-L^" 2 (-Pl)n(-7-Pl + l) w ( ir+ l, 

(A.2) 



r [7 + pi] r [-P2] \V^J 

m ( t 



r [7 + in] r [-pi] 

whereas for the far region, equation (5.22), we get 



F , t s = ~ r [pi - P2] £[g rr-, t yi + P2) n {u + p2) r 



T[a + pi]T[p + pi] vv ' (7i-l) 



(v^<) pl T[s _ a _ pi] r [s-p-kY (a - 3) 

The comparison between the two equations leads us to the conclusion that P2 — P\ = —n 
doesn't constitute a solution since the matching of the two solutions is not possible. Finally, 
we should recall that we have assumed neither a nor b of the problematic Hyper geometric 
functions in equations (5.24) and (5.22) are equal to a negative integer smaller in absolute 
value then n. We now consider the case where this assumption does not hold. 

In this case, the Hyper geometric functions in question no longer diverge and we can 
use the previous expansion F[a,b,d;s] ~ 1 + O(s). However, comparing equations (5.22) 
and (5.24) with (A.2) and (A. 3), respectively, we see that if a or b fulfill this requirement 
only in one side of the equation, i.e only in the far or near region, we will again be unable 
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to do the asymptotic matching since the powers of t on both sides of the matching equation 
are different. Now we can re-examine equation (5.26) in the case a = —m or b = —m and 
n < m for both the far and the near region. For the far region we have 



while for the near region 



a far =/3 + P2 

b far =a + p 2 (A. 4) 

ir = — Pi 

,r = -pi-7 + l- (A.5) 



Now we can rewrite equation (5.26) as 

2 



(v[n]) 2 v[ 7 + p 2 ]r{ ] V [a far ] T [b 



T[a + Pl ] T[f3 + Pl ] (T [-„]) 2 T [7 + pi] T [-p 2 ] ( ^ (A ' 6) 

We observe that, for this equality to be fulfilled, we need to have at least one more divergent 
gamma function in the denominator then in the numerator of the LHS since the RHS goes 
to zero. We will show this doesn't happen, assuming throughout our considerations the only 
divergent term in the denominator is (r [—n]) 2 since the other cases were already treated. 
Thus, we need to check only the cases where there is only one divergent gamma function 
in the numerator. This can only happen when aj ar , bf ar are equal to —mi, —m 2 > — n and 
bnear also fulfills this requirement but a near doesn't fulfill it. 

We must assume then that (l)pi — pi = n and (2) — — p\ > —n. First, we observe 
the equality pi = —pi — 7 is true for the specific parameters of our equation. This implies 
that pi + 1 = — m and hence pi G Z and from (1) also pi £ Z. Next, inserting 7 = 2p + 1 
into assumption (2), we find —pi > —n + 2p. However, we also need to have pi < —n, 
since we showed it is an integer. It follows that 2p < 0, and since pi G Z we have 2p G Z. 

Using (1) we shall now prove that the in this case Tpy + pi] is divergent. 

p 2 + -i = l + n-pi + 2p<l + 2p<l 

since both 7 and pi are integers, obviously 7 + P2 is an integer, at most equal to zero. This 
ends our proof. 

B Analytic expressions for the radial part of the eigenfunction 

First we present the explicit expression for the radial part as a function of the separation 
constant A in the region near infinity , corresponding to ^fzi < r < 1, 

R(r) = (-1 + r f r K~ 1+ \fi~^R^ 2 Fx [A,B,3,l-r] (B.l) 

-1 + a 2 ) Jl + a 2 + %a 2 m^ 




aX 

a + a(4 + aa) (2 + aa)(a + a (4 + aa)) 



(B.2) 



aX (-l + a 2 )^l + a2 + ^aV» 

H : „„\r„ i „ia i „„u — • ( B - 3 ) 



a + a(4 + aa) (2 + aa)(a + a(4 + aa)) 
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where we have used our previous definition of 



(B.4) 



x 2 + (1 + a) 2 ' 

For the case of the first stable mode of a black hole with a = 2 we have obtained A 
analytically and can thus simplify the above expression to the form 



R(r) 



2(1 + a)r 5" ( (1 + o)ri+° + r i+° (-1 - a + (-1 + o)m, 
(-1 + a) (-1 - a + (-1 + a)m^) 



+ 



2(1 + a)r 1+ 3P+=T 
1 + a + (1 - a)m^ 



(-l+a)« 



(B.5) 



In the near horizon region , corresponding to < r < y/z2, we shall not present here 
the general expression but rather only that for the case of a = 2 where we obtained A 
analytically. It is given by 



R(r) = 



. \-2n -p T-i (-l+a)m ] _ [ - (-l+o)mjl 2 o- [(-l+o)irmil n 

2 ( 1 + a ) r [ 1 i+7- £ J r [- 1 -"+ — I+^J Sm [ ThT^J z 2 



1 - 2n + 



-l + a)m^ "I 
1+a J 



2^1 



— n, — 1 — n 



(-1 + aW x _ 
1 + a 



— l + a)m^ 

Z2 



1 + a 



(r-l)V 1_n+< At+iT*x 
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